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1 Introduction 

I describe in this talk recent work on bosonization developed in colaboration 
with my friends Nino Bralic from Universidad Catolica de Santiago de 
Chile, Cesar Fosco from Centro Atomico de Bariloche, Eduardo Fradkin 
from the University of Illinois at Urbana Champaign, Jean Claude Le 
Guillou from LAPP-Annecy and ENS de Lyon, Enrique Moreno from 
City Colege and Baruch College of the CUNY, and Virginia Manias and 
Carlos Nunez from the University of La Plata. More details can be found 
in references 0-0. 

Bosonization is a mapping of quantum field theories of fermion fields 
onto equivalent theories of boson fields. Well-established in 1 + 1 dimensions, 
bosonization constitutes there one of the main tools available for the study 
of the non-perturbative behavior of both quantum field theories || and of 
condensed matter systems 0. In dimensions other than 1 + 1, much less 
is known. In this talk I will precisely discuss the issue of d > 2 bosoniza- 
tion within the path-integral approach making special emphasis in the d = 3 
case, of particular interest in condensed matter problems. The bosoniza- 
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tion approach that I will present corresponds to the path-integral framework 
and basically establishes an identity between the generating functional of 
the fermionic theory and the generating functional of the equivalent bosonic 
theory. From this relation, the recipe for bosonization of fermion currents 
is derived and the current commutator algebra is presented. To complete 
the bosonization program one should also calculate the energy-momentum 
tensor algebra. This is not done in this work. 

Related and unrelated approaches to d > 2 bosonization can be found in 
Kj-[g3 . 



2 The method 

The method is straightforward. We shall basically consider the case of free 
fermions but we shall also discuss an interacting (Thirring) model. 

One starts from the fermion Lagrangian for N massive free fermions in d 
dimensions, 

L = + m)ij) (1) 
The corresponding generating functional reads 

Z fer [s)=[vipV'ipexp[-[ d d xi>(i$ + ^ + m)ip] (2) 

J JM 

where s M is the source for fermion currents, s M = s a ^t a with t a the generators 
of a group G and M a <i-dimensional manifold. 

Our derivation heavily relies on the invariance of the fermion measure 
under local gauge transformations h(x) G G with G the group of continuous 
maps M —>■ G. This ensures that 

Z[s h ] = Z[s] (3) 

with 

si = h'h^h + h^d^h (4) 
Evidently, fermions can be integated in (|J), 

Zf er [s] = det(i$ + $ + m) (5) 
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and the determinant in the r.h.s. of eq.(|5]) will be used to introduce an 
auxiliary field taking values in the Lie algebra of G through the trivial 
formula 

Z fer [s] = J Vbp8V>(bp - Sft ) det(0 + y + m) (6) 
Now, it will be convenient to replace the delta function in (]£]) as follows 

Here we have used that the equation 

Mb]=fM (8) 
has for s M ^ the unique solution 

K = ^ (9) 
and A [b] is a Faddeev-Popov-like jacobian, 

A[b] = \det(2e fllfi2 ... fld D fll [b})\ (10) 

with D^b] the covariant derivative, 

D,[b] = d,+ [b„ ] (11) 

We do not consider for the moment Gribov like problems that could arise for 
certain manifolds and groups. Concerning the delta function in the r.h.s. of 
eq.flT]), n depends on the space-time dimensions according to n — d(d — l)/2 
since one needs for enforcing eq. @ as many 5-functions as independent com- 
ponents the curvature has. 

It is at this point that the bosonic field whose dynamics will be equivalent 
to that of the original Fermi field comes into play. We introduce it as a 
Lagrange multiplier A^ 3 „. M enforcing the ^-function in the path-integral ffl), 



r /er [s] = jvb^VA^ d det(i$ + y + m)A[b] 



cxp 



\trjd xA^ 3 ^^ d e fllfl2 ,,, fld (f^ lfl2 [b] fni(i 2 [ s )) 



;i2) 
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Here A is a constant which can be adjusted so as to obtain an adequate 
normalization for the currents. Now, we rewrite eg. (12) in the form 

^fer ■■^d eyL P(Sbos[A]) exp ( ^trj d d xA^ 3 _^ d e fll ^ 2 _^ d f^ lfl2 [s]\ 

(13) 

where the bosonic action is defined as 

exp (S bos [A] ) = J Vb„ det + m)A [b] 

exp (\tr J d d xA / j t3 _ M ^ d £ fllfM2 _ tJid 

(14) 



Formulae (|13])-(|T4]) constitute our basic bosonization recipe: eg. (p"3D allows 
to compute fermion current correlation functions in terms of the bosonic field 
A and eg. ([14]) gives the bosonic action defining the dynamics of A. It can 
be now appreciated in what sense we consider our bosonization recipe exact: 
we have arrived with no approximation to a bosonization recipe of the form 

^/^2Ae wla ... w 9 M2 AJ 1 ... w (15) 

However, except in d = 2 dimensions where we know how to compute exactly 
the fermion determinant appearing in (|T4]) and to resolve the path-integral 
defining Sb os , one should appeal to some approximation scheme to evaluate 
the bosonic action accompanying recipe flT5|). This means that only in d = 2 
dimensions the complete bosonization recipe is exact. 

It should be stressed that the bosonization recipe fll5| ) should be taken 
as illustrative of the bosonization since the rigorous eguivalence between the 
fermionic and the bosonic theory is at the level of the generating functional 
Zf er [s] of Green functions. It is from Zf er [s] written in the form ([U|) that 
one has to compute current correlation functions in the bosonic language. 
Note also that in writing recipe (|T5D we have ignored terms which are non 
linear in the source. Although correlation functions of currents acguire a 
contribution from this terms, this contribution is irrelevant in the calculation 
of conmutator algebra since they have local support. This can be easily seen, 
for example, using the BJL method (see |5| for a discussion and JF| for the 
application of the BJL method within the present bosonization approach). 
Concerning the bosonic field A" „ note that it corresponds to scalar fields 
in d = 2 dimensions (see ref.0 for details on how to make contact with the 
usual bosonization rules), to a vector field in d = 3 dimensions and to an 
antisymmetric (Kalb-Rammond) field in d > 3 dimensions 13], pf. 
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3 The abelian case in d = 3 



The Abelian case in 3 dimensions is particularly simple. To begin with, the 
bosonization recipe for the fermion current reads 

_ i 

V'TmV' -> ±T-Cp,ad u A a (16) 

where we have chosen A so as to make contact with the normalization em- 
ployed in ref.@. Concerning the bosonic action, A [b] is trivial so that eq.flljj) 
simply reads 

exp(S bos [A)) = J Vb^ det(i$ + h 1 + m) exp (tj^ tr J d z xA^e^ a f va [b] 
or, calling 



(17) 

-logdet(i# + J(+m) = / d 3 xL[b] (18) 



we can write 



exp(S bos [A]) = J Vb„exp(-S e ff[b,A]) (19) 
where 5 e // is defined as 

S eff [b, A} = J d 3 x(L[b] ±^tr A^af^b]) (20) 

(The double sign in eqs.([T7| ) - (|20|) is included for convenience, see the discus- 
sion below) 

Being in general L[b] non-quadratic in b one cannot path- integrate in 
flT7D so as to obtain 5& os [A]. We shall see however that there is a change of 
variables allowing to decouple A^ from b^ in S e ff[b, A] so that one can control 
the A^ dependence of S& os [A] without necessity of explicitly integrating over 
b^. Let us define a new variable b'^ through the equation 

b ll = (l-0% + OA ll + V ll [A] (21) 

where V^A] is a gauge invariant function of A^ so that b'^, the variable which 
will replace b^, transforms as a a gauge field. 9 is an arbitrary parameter to 
be adjusted later. The idea is to choose so as to decouple b'^ from A. This 
amounts to impose the following condition 

62Seff - (22) 



6VJx)6A v (y) 
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which in terms of V p reads 



2i\e paa d"5(x — y) + 9 y J d 3 z 



5b p (y)5b a (x) 



J \J 5b p (u)5b a (x) J 5A p (y) 

To go on we need an explicit and necessarily approximate expression for L [b] . 
If we are interested in the large-distance regime of the bosonic theory we can 



use the result first [p4j| -[p5 for the d = 3 fermion determinant as an expansion 



in inverse powers of the fermion mass 



L[b) = T ^e pua b p d u b a + ^-JlM + 0{\) (24) 
Sit 24Tr\m\ p m z 



The first term in (Ell) is the well-honored Chern-Simons action introduced in 
7M as a way of generating a mass for gauge fields in three dimensions. The 



double sign in this term is originated in a regularisation ambiguity charac- 
teristic of odd-dimensions (see ref.|2(J). The second one corresponds to the 
leading parity-even contribution to the fermion determinant. One easily sees 
that if one tries for V^A] the functional form 

V ll [A}=i^-e pua f ua [A] (25) 



one gets, from the decoupling equation (|23"D, 

C = ±1/3 (26) 

Then, if for simplicity one chooses 9 = — 1, the bosonic action for A p can be 
easily found to be 

S bos [A] = ±i- / d\e pua A p d v A a + — \— f d\f\A] + 0(l/m 2 ) (27) 

o7T J z47r|m| J ^ 

One can in principle determine, following the same procedure, the following 
terms in the 1/m expansion of Sb os by including the corresponding terms 
in the fermion determinant expansion. This result extends that originally 
presented in ref. [|TJ. It shows that the bosonic counterpart of the three di- 
mensional free fermionic theory is, to order 1/m 2 , a Maxwell-Chern-Simons 
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theory which is equivalent, as it is well-known, in turn to a self-dual system 



Alternatively to the 1/m determinant expansion, one can consider an 
expansion in powers of 6 M retaining up to quadratic terms. The result can 



be written in the form ||17 
L[b\ 



i 1 
-j£^ab^Pd v b a + — — f^[b}Qf pu [b] 



(28) 



where P and Q are functionals to be calculated within a loop expansion, 



P = P 



2 



Q = Q 



(J 2 



(29) 



Details of the calculations of P and Q and results within the loop-expansion 
can be found in refs. |T7|,P^|. 

In order to decouple the 6 M field one again proposes a change of variables 
like in fl2"T| ) but now trying for V p the (gauge-invariant) functional form 



with 



% % 
= — £uuaPfva[A} = 2 — e, lua Rd u A a 
m m 



d 2 



R = R 

One finds, from the decoupling conditions 

S 2 S eff [b',A] 
5A p (y)5bUx) 



(30) 



(31) 



d 2 



(1-9) te paa [2X + 6P- 2—QR d Q 5(x - y) 



77T 



+ - (\0Q - PR) (d p d a - 5 p(J d 2 ) S(x - y) 



m \2 







9 being here a functional of d 2 /m 2 . The solution of this equation is 

Q P 



R = -A 



(P 2 - ^ 



9 = -2A 



P 2 - Kq 2 



(32) 



(33) 
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With this choice, the change of variables decouples the b^ integration so that 
one can finally get the bosonic action for which now reads 



S bos [A] = Jd 3 x i-{2\f l -e^ a A, 



P 



+(2A) 2 — MA] 



(p 2 - 
Q 



■d u A a 



(pi _ 



(34) 



This result coincides with that found in ref.||17||, obtained by a direct func- 
tional integration on b^. As it was proven in this last work, it corresponds 
for massless fermions to the bosonization action proposed in ref. |15| since in 
the m — ► limit eq.(|34]) takes the form 



2 r 3 1 1 % 

Sbos — —id x{-Fn V —j==Fau — -tau\Aad u A) 

7r J 4 v —a 2 2 



(35) 



4 Interacting models 

One can apply the bosonization approach described above to analyse inter- 
acting fermionic models. Let us consider for example the Thirring model 
with a current-current interaction Lagrangian of the form 



9 Z 



J int 



2N 



where ip % are N two-component Dirac spinors and the U(l) current, 



(36) 



(37) 



The coupling constant g 2 has dimensions of inverse mass. (Although non- 
renormalizable by power counting, four fermion interaction models in 2 + 1 
dimensions are known to be renormalizable in the 1/N expansion |30| . ) One 
can directly apply the bosonization recipe found in the precedent section to 
this interaction Lagrangian. Making a choice of A so as to coincide with the 
normalization in U, this meaning 



'^e dA 



(38) 
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one has 

L int - £-F% (39) 

Alternatively, one can eliminate the quartic fermionic interaction by intro- 
ducing an auxiliary field via the identity 

exp( J £jj"jA ) = / Va » ex P[- / (\ af1 ^ + -jfifajfo ] (40) 



N 

and then proceed to integrate fermions as in the free case thus obtaining a 
determinant in which the a M field can be eliminated by a shift — > — a M . 
One confirms in this way that the bosonization recipe ( |39"D is correct so that 
the three dimensional Thirring model is equivalent, in the 1/m approximation 
to a Maxwell-Chern-Simons model. To leading order in 1/m we can then 
write (after rescaling the field A^) 

1 27T 

Lth - -Flu ± i^e^A^ a A v (41) 

We can give now a first application of the bosonization formulas and, in 
this way, explore their physical content. The Lagrangian in ([|l]) has a 
Chern-Simons term which controls its long distance behavior. It is well 
known [pi], j3l"H that the Chern-Simons gauge theory is a theory of knot in- 
variants which realizes the representations of the Braid group. These knot 
invariants are given by expectation values of Wilson loops in the Chern- 
Simons gauge theory. In this way, it is found that the expectation values of 
the Wilson loop operators imply the existence of excitations with fractional 
statistics. Thus, it is natural to seek the fermionic analogue of the Wilson 
loop operator Wr which, in the Maxwell-Chern-Simons theory is given by 



W r = (exp{z^— f A^dx^}) (42) 
9 J r 



where V is the union of a an arbitrary set of closed curves (loops) in three 
dimensional euclidean space. Given a closed loop (or union of closed loops ) 
T, it is always possible to define a set of open surfaces £ whose boundary is 
T, i.e. T = <9X. Stokes' theorem implies that 



(exp{z— I A^dx*}) = (exp{z^ f dS^ x d u A x }) 
g Jr g 
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(exp{z— f d 3 x e^ x d u A x b x }) (43) 



<N 
9 

is an identity. Here b\(x) is the vector field 

b\(x) = n x (x)5x(x) (44) 

where n\ is a field of unit vectors normal to the surface £ and Ss(x) is a delta 
function with support on E. Using eq.(|3"8|) we find that this expectation value 



becomes, in the Thirring Model, equivalent to 

W v = (exp{z — / dx^}) MCS = (exp{ / dS^^}) Th (45) 



'N 
9 

More generally we find that the Thirring operator We 



W s = (exp{g ^ dS^Y^hh (46) 

obeys the identity 

r - Vn r 

(exp{g / dS^Y^Th = (exp{iq & A^}) M CS (47) 

for an arbitrary fermionic charge q. 



The identity (47) relates the flux of the fermionic current through an 
open surface £ with the Wilson loop operator associated with the boundary 
T of the surface. The Wilson loop operator can be trivially calculated in the 
Maxwell-Chern-Simons theory. For very large and smooth loops the behavior 
of the Wilson loop operators is dominated by the Chern-Simons term of the 
action. The result is a topological invariant which depends only on the linking 



number vy of the set of curves T |Tl], [J!]]. By an explict calculation one finds 

f - Nq 2 
(exp{g / dS^^tp^Th = exp{=Fii^-r— } (48) 

JT, 87T 

This result implies that the non-local Thirring loop operator We exhibits 
fractional statistics with a statistical angle S = Nq 2 /8n. The topological 
significance of this result bears close resemblance with the bosonization iden- 
tity in 1 + 1 dimensions between the circulation of the fermionic current on 
a closed curve and the topological charge (or instanton number) enclosed in 
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the interior of the curve [p2 |. From the point of view of the Thirring model, 
this is a most surprising result which reveals the power of the bosonization 
identities. To the best of our knowledge, this is the first example of a purely 
fermionic operator, albeit non-local, which is directly related to a topological 
invariant. 



5 Current algebra 



We have seen in precedent sections that in 3 dimensional space-time the 
fermion action bosonizes, in the large m limit, to a Maxwell-Chern-Simons 
theory. Now, the gauge invariant algebra of such theory has been studied in 
refs. , ||28|| . One has for instance with our conventions, 



[Ei(x,t),B{y,t)] = -3\m\e l3 d^ 2 \x~y) (49) 

If one now relates the electric field E^ = F iQ and the magnetic field B = 
6ijdiAj to the fermionic currents through the bosonization recipe for the 
fermion current, 

Jo - 4=B (50) 

V47T 



'An 

then, the resulting fermion current commutator algebra is not the one to 
be expected for three-dimensional free fermions. Indeed, the d = 3 fermion 
current algebra should contain an infinite Schwinger term |33]-|]35| which is 



absent in eq. (|49"D . The point is that calculations leading to a bosonic theory 
of the Maxwell-Chern-Simons type are valid only for large fermion mass 
while calculation of equal-time current commutators imply, as we shall see, 
a limiting procedure which cannot be naively followed for large masses. 

Since the exact bosonic partition function is much too complicated to 
handle, a possible strategy is to use the quadratic (in auxiliary fields) ap- 
proximation mentioned in the precedent section working with an arbitrary 
(not necessarily large) mass so as to obtain a bosonized version of the origi- 
nal fermionic model in which the equal-time limit can be safely taken. One 
should then compute current commutators for this bosonized theory, and 
test whether they coincide with those satisfied by fermionic currents in the 
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original model. Details of this calculation can be found in M. I will just 
sketch here the principal steps leading to the consistent equal-time current 
commutators in the bosonic language. 

As explained above, within the quadratic (in 6 M ) approximation, one can 
write the fermionic partition function in terms of the bosonic fields in the 
form fll? | 

Zfer = J DA^exp - j d z x(-F,jvCiFp V -^ApCzepvxdvAx 

+is fJ ,e flv xd v Ax)] (52) 

with C\ and C2 now given through their momentum-space representation C\ 
and C2 

1 F(k) 



C^k) 
C 2 {k) 



An k 2 F 2 (k) + G 2 (k) 

1 G{k) 
4tt k 2 F 2 {k) + G 2 {k) 



(53) 
(54) 



and F{k) and G(k) given by |HJ 



F(k) 



I m I 
Ank 2 



1 - 



k 2 
y Am 2 ' 



Am 2 ■ n , 4m 
- arcsmll + 



A- 2 



G(fc) = — 



Hi 



Q 



An 2n\k 



arcsin(l + 



Am 



Being quadratic in A^, eq.([52|) can be integrated leading to 



in 



d 3 xd 3 yd u s l j i (x)e l j iU xF>xp(x,y)d a s T (y)e 



par 



J. 

Zbos[s] = [detD^] 2 exp 

(55) 

where D 1 is just the propagator of the bosonic action which, in the Lorentz 



gauge we adopt from here on, reads 
d 3 k 



(2tt)* 



[P(k)g,j, v + Q(k)k^k v + R{k)e^ va k a ] expik(x - y) 



(56) 
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with 

pw= ^5+W) = ^ (fc) (57) 

o( k ) = dl{k) (A^L) 2 (58) 

R{k) = _f 2(fc) _„ (59) 

Let us briefly recall how one can compute current commutators within the 
path-integral scheme using the so-called BJL method fl36|-|38|. To this end 
we define the correlator 



G„(*,y)- PlosZ,M 



8s^(x)8s v (y) 



(60) 



s=0 



from which one can easily derive equal time current commutators using the 
relation 

< [jo(x, t),ji(y, t)] >= lirn [G w (x, t + e;y,t)- G 0i (x, t-e;y, t)] (61) 

The current commutator evaluated using eqs.(]5"0"|)-(|6"T|) corresponds to Zf er [s] 
written in terms of bosonic fields. That is, eq.(|6T|) gives the equal-time 
commutator for the bosonic currents = (l/\f4xr)e fJ , l , a d v A a . This result 
should then be compared with that arising in the original 3-dimensional 
fermionic model for which = —iip'y fl ip [|35j . 

Starting from eqs.(p5|)-(|56]) and using the BJL method we get, after some 
calculations, 

G^fay) = -—e mp 6 vl3u d a dpD^ (62) 

or 

1 f d 3 k 

G^fay) = — J j^-^[P(k)(k 2 g^ - k^k v ) + ^Ri^e^k^ exp[ik(x - y)] 

(63) 

With this, we can rewrite eq. (|6~l"D in the form 

< \jo{x, t),ji(y, t}>= lim F(x - y) (64) 
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with 

r(x) = -2i 
where we have written (/c M ) 



d 3 k 
(2^)3 
(fc , k. 



kokiSm(k e)F(k) exp (ik.x) (65) 
, i — 1, 2. It will be convenient to define 



(66) 



In terms of this new variable and using the explicit form for F(k) given in 
0)0 > w ith A; = (ko + k 2 ) , integral V becomes 
1 1 „ r d 2 k 



P(x) 



^Tx 2 \m\ e 2 



0; 



expi(k.x) / dk' Q k' sin k' f(y) 

(27TJ Z JO 



where 



M 

and we have defined 



1 

y 



arcsm • 



2/ 



"~ y/l + (1/y) 
k 2 k' 2 + e 2 k 2 



4e 2 m 2 



(67) 



(68) 



(69) 



4m 2 

One can now see that y — > oo for e — > and fixed m. Then, expanding in 
powers oil/y one has f(y) ~ Tv/(2y/y) and then using distribution theory to 
define the integral over /c one finds 



— lim -di5^ (x — y) 



(70) 



This result for the equal-time current commutator, evaluated within the 
bosonized theory, shows exactly the infinite Schwinger term that is found, 
using the BJL method, for free fermions in d = 3 dimensions ||35|| . As it 



happens in d = 4 dimensions [34], we see from eq.(70) that the commutator 
at unequal times is well defined: divergencies appear only when one takes the 
equal-time limit. 

One can evaluate also the next order vanishing in the equal-time current 
commutator so as to compare it with the result from the original fermion 



model reported in the literature |p5| . The answer is J3| 
<\jo{x,t),ji{y,t))> = -^-\im(-di5 (2 \x-y) 

-UAm 2 di5 {2 \x - y) - -diA6 {2) ^ 



x 



(71) 



15 



where we have defined 

1 f°° 1 

A = io *°^fin* (72) 

In order to compare with ref. KJ] where current commutators were computed 
using dimensional regularization, we define, coming back to the original vari- 
able k = k' /e 

1 f 1 

A[d) = - d d ^ 2 k -^sink e (73) 

so that A[d = 3] = e/A. One can now perform the analytically continued 
integral to find, near d = 3, the behavior 

A[d] ~ -e x ill (74) 



The same ambiguous result for free fermions is obtained in ref.|35[ near d = 3. 
This ambiguity can be however removed, the pole in dimensional regulariza- 
tion corresponding as usual to a logarithmic divergence. It is also interesting 
to note that if one uses the nice approximation F appr for F proposed in 



ref. [17], one can well check the correctness of our previous analysis 



From the analysis above, we see that not only the infinite Schwinger term, 



analogous to that arising in d = 4 |34| is obtained in the bosonized version of 



our d = 3 fermion theory but also the mass-dependent second term as well 
as the triple derivative third term, both vanishing in the equal time limit. 
Our analysis should be compared with that in ref.|39[] where the fermionic 
commutator algebra is inferred from the Maxwell-Chern-Simons algebra for 
electric and magnetic fields using a bosonization recipe which is valid in the 
large mass limit. One can see that in the large mass regime, terms depending 
on the product em = A will produce ambiguities according to the way both 
limits (e — > and m — > oo) are taken into account, a problem which is not 
present in the limit of small masses. To see this in more detail, let us come 
back to (|65|) and consider the case in which A is kept fixed while e — > (so 
that m — > oo). In this case, taking the limit before integrating out k' , one 
finds for F 

F(x) ~ \m\h(\)di5 {2 \x) (75) 

where 

h(\) = — dzzsm(2\z)f{z) (76) 

27T JO 
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with / given by eq.fl68|). Let us note that using the approximate F appr of 



ref. ||17|| and taking the limit after the exact integration over k' , we recover the 
same behavior Q75|). We see that for A = em fixed, h just gives a normalization 
factor so that one reproduces from P in the form (|75|) a commutator algebra 
at equal times and large mass that coincides with that to be infered from a 
Maxwell-Chern-Simons theory, 

< [jo(x,t),ji(y,t} > — ► c|m|^5 (2) (f - y) (m -> oo) (77) 

with c a normalization constant. Again, currents appearing in the l.h.s. of 



eq. (|77|) are bosonic currents which can be written in terms of the electric and 



magnetic fields thus reproducing the MCS gauge invariant algebra |24fl.[p£ 
One should note however that the free fermion - MCS mapping is valid in 
the large mass limit of the original fermionic theory, this meaning the large- 
distances regime for fermion fields. Since current commutators test the short- 
distance regime, one should not take the MCS gauge-invariant algebra as a 
starting point to reproduce the fermion current commutators. 
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6 The non-Abelian case in d = 3 



In three dimensional space-time, the bosonic action ([14]) takes the form 
exp(-S bos [A}) = J Vb^Vc^Vc^exp (-tr Jd 3 x (L[b\ ± -^e^ va c^D u [b}c a 



Here ghost fields and were introduced to represent the Faddeev-Popov 
like determinant A [b\. Again, we have written 

tr Jd 3 x L[b] = - logdet(z# + m + f) (79) 



and we have chosen the arbitrary constant A appearing in (|14l) so as to make 
contact with the conventions of ref.0, A = jj-. Moreover, we have shifted 
the bosonic field — > — b^ (this amounting to a trivial Jacobian) for 
reasons that will become clear below. 

It was observed in ref. || that when L[b] is approximated by its first term 
in the 1/m expansion, a set of BRST transformations can be defined so 
that the corresponding BRST invariance allows to obtain the (approximate) 
bosonic action. We shall explicitly prove here that this invariance is present 
in ( |78|) where no approximation for L[b] is assumed. To this end, we introduce 
a set of auxiliary fields (taking values in the Lie algebra of G) , / and x so 
that one can rewrite (|78|) in the form 

exp(-S bos [A]) = JvbvVc»Vc»VhvVlVxexp{-Seff[A,b,h,l,c,c,x)) (80) 
with 

S ef f[A,b,h,l,c,c,x] = trJd 3 x(L[b-h)±-^e flva c fl D u [b]c a 

T^((A, - bJJJb] + lh\ - 2xKc,)) (81) 



where x is an anti-ghost field. Written in the form (plD , the bosonic ac- 
tion has a BRST invariance under the following nilpotent off-shell BRST 
transformations 

5c M = Af, - 6 M , 8b ^ = c M , 8A^ = c M , 5c M = 0, 8x = I 
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5h„ = c„, 51 = (82) 

In view of this BRST invariance, one could add to S e ff a BRST exact 
form without changing the dynamics defined by S6 OS [A|. Exploiting this, we 
shall see that one can factor out the dependence in the r.h.s. of eq.(|81"|) 
so that it completely decouples from the path-integral over auxiliary and 
ghost fields exactly as we did in the Abelian case. Although complicated, 
this integral then becomes irrelevant for the definition of the bosonic action 
for A^. Indeed, let us add to S e ff the BRST exact form SG, 

S eff [A, b, h, I, c, c, x) -> Seff [A, b, h, I, c, c, x] + SG[A, 6, h, c] (83) 

with 

G[A, b, h, c] = T^-tr d\e^ a c^H ua [A, b, h] (84) 

107T J 

and H ua [A,b,h] a functional to be determined in order to produce the de- 
coupling. Then, consider the change of variables (analogous to ( pT|) for the 
Abelian case) 

b ll = 2b'^-A lx + V lx [A] (85) 

where V^A] is some functional of A^ changing covariantly under gauge trans- 
formations, 

V»[A°} = g- l V,[A]g (86) 

so that b'^ is, like A^ and b^, a gauge field. Integrating over I in ( |8"UD and 
imposing the resulting constraint, = 0, one sees that if one imposes on 
H ua [A, b, h] the condition 



(87) 



6M(y) 6A a (y) 5h a (y) 



then, when written in terms of the new b'^ variable, the ghost term becomes 
S g host[b', c, c] = ±^-tr (fxs^aC^D^b'jca (88) 

o7T J 

so that its contribution is still A^ independent. Then, we can write the 
effective action in the form 

S e f f [b',A]+S 9host [b' 7 c,c] (89) 
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with 

S eff [b',A] = S[b,A] 

= trjd 3 x (L[b]T^(A^-b^(%[b]+*H^[A,b,0]^ (90) 



where = e^ a H va . 



Condition ( p7|) made the ghost term independent of the bosonic field A^. 
We shall now impose a second constraint in order to completely decouple the 
auxiliary field b'^ from A^ in S e ff. Indeed, consider the conditions 



5A*(y)5V>(x) 

In terms of the original auxiliary field 6 M these equations read 

5 2 S[b,A] 8*S[b,A] | r, u 6*S[b,A] 6Vj(u) _ Q (g2) 



5A;(y)5b» a (x) 6b a p (y)6b b a (x) J 5b%(u)5b b a (x) 5A«(y) 



Eqs.(^2j) can be easily written in terms of L, H and V as a lengthy equation 
that we shall omit here. 

The strategy is now as follows: once a given approximate expression for 
the fermion determinant is considered, one should solve eq.(]9"2]) in order to 
determine functionals V in eq.(^) and G in eq.(^), taking also in account 
the condition (|87|) . In particular, if one considers the 1/m expansion for the 
fermion determinant, equations fl87|) and (|92|) should determine the form of 
V and G as a power expansion in 1/m. In ref.|25j the 1/m expansion for the 
fermion determinant was shown to give 

lndet^ + m + Jfi = ±-^-S C s[b} + I PC [b] + 0(d 2 /m 2 ), (93) 

where the Chern-Simons action Scs is given by 

f 2 

Scs[b] = e^x tr J d 3 x (f^b x - -b^bx). (94) 

Concerning the parity conserving contributions, one has 

Ipc[b] = tr j d\ rf» v + • • • , (95) 
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To order zero in this expansion, solution of eqs.(|S7D,(|9"2"D is very simple. 
Indeed, in this case the fermion determinant coincides with the CS action 
and one can easily see that the solution is given by 

VW[A) = (96) 
G® [A, b, h, c] = ±JLtr Jd\ c, (i %[A] + ~ %[b] - TD, a [A}h a ) (97) 



With this, the change of variables ( pop takes the simple form 

b, = 26; - A, (98) 
and the decoupled effective action reads 



Si%[b, A, c, c] = T --(2Scs[b'} - S CS [A\) + S ghost [b'\ (99) 

±07T 

We then see that the path- integral defining the bosonic action Sb os [A], factors 
out so that one ends with a bosonic action in the form 

SfM] = ±^-S cs [A} (100) 

as advanced in 0,0. Let us remark that in finding the solution for G one 
starts by writing the most general form compatible with its dimensions, 

G (0) [A, b,h,c) = tr J d 3 x e^aCp (dib u A a + d 2 A u b a + d 3 b u b a + d±A v A a 

+d 5 b„h a + d G h„b a + d 7 A u h a + d^h u A a + d 9 d v A a + diod u b a + dnd v h a ) 

(101) 

All the arbitrary parameters di are determined by imposing the conditions 
fl57D and Q9~2|) with *H fl transforming covariantly (as does) under gauge 
transformations which leads, together with a gauge invariant action, to the 
solution (P7|) . 

To go further in the 1/m expansion one uses the next to the leading 
order in the fermion determinant as given in eq. (|9"BD . Again, starting from 
the general form of G and after quite lengthy calculations that we shall not 
reproduce here, one can find a unique solution for V M and H ua leading to a 
gauge invariant action, 

^ (1) ^]=±^%[^] (102) 
21 
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[A, b, h, c] = G (0) [A, b, h, c] =F gg-^^ Jd 3 x c p e^ a e„ pa 
( \ [ - h] + 3f pa [b -h}- 2D p [A - h](A a - b a ) , (A* - b a ) } 
+ 4[f pa [A~h], K) ) (103) 
The corresponding change of variables fl55|) takes now the form 

& M = 2&; t -A M ±£-%[^] (104) 
and the decoupled effective action reads 

S&fc A, c, c] = Si%[b, A, c, c}+tr Jd\ (g^/iM + 24^&H) (105) 

so that one can again integrate out the completely decoupled ghosts and b' 
fields ending with the bosonic action 

Sll][A] = ±^S CS [A] + ^-tr Jd*x fa[A) (106) 

This result extends to order 1/m the bosonization recipe presented in refs. 
1,1- 

In this way, from the knowledge of the 1/m expansion of the fermion 
determinant one can systematically find order by order the decoupling change 
of variables and construct the corresponding action for the bosonic field A p . 
One finds for the change of variables 

b p = 26; - A, ± ^%[A\ + ^C^D P [A\U P [A\ + ... (107) 

Here is a (dimensionless) constant to be determined from the 1/m 2 
term in the fermion determinant expansion, which should be proportional to 
*f ll D p f pp . Evidently, finding the BRST exact form becomes more and more 
involved and so is the form of the bosonic action which however, can be 
compactly written as 

S bos [A}=trJd 3 x (l[-A + V[A}} 

T-^(2A P - V,[A})(J p [-A + V[A\] + *H»[-A + V[A], A, 0]))(108) 
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Let us end this section by writing the bosonization recipe for the fermion 
current accompanying this result for the bosonic action. From eq.(|T5D we 
have, in d = 3 

#7^' - ±^ va d v A a a (109) 

7 Wilson loops 

I will transcribe in this section some results obtained in || concerning the 
evaluation of Wilson loops in the framework of our bosonization approach. 
In the Chern-Simons theory, they measure topological invariants determined 



by the linkings of the loops and by the topology of the base manifold |Jl 
For one loop T, 

W[F] = ti P exp(i j> da? (110) 

where P denotes the path ordering of the exponential, and the trace is taken 
in the representation carried by the loop. According to the bosonization 
prescription, to relate this operator to the fermionic theory we must express 
W\T] in terms of the field strength rather than the potential A^. In 
the abelian case this can always be done by means of Stokes theorem. As 
discussed in ref . || , this leads to an explicit mapping between abelian Wilson 
loops and non-local fermionic operators. Hence, in this way, the latter are 
related to the linking of loops and thus probe the generalized statistics of the 
external particles that propagate along those loops. One way to extend that 
analysis to the non-abelian case is to use the non-abelian extension of Stokes 
theorem developed in [flOfl . For an arbitrary loop V = <9£, the boundary of a 
surface S, one has 

W[dZ]=tTP t exp{ij o dtj ds ^ - H/- 1 [ s S(t) ]F^(S(t,g))H/[ s S(t)o]} 

" ( ' S (111) 
Here £ is looked upon as a sheet, that is, a one parameter family of paths 
parametrized by t, < t < 1. For each t, T,(t) is a path, itself parametrized 
by s, < s < 1, with fixed end-points: dY>(t,s)/dt = at s = 0, 1. For a 
given t, s T,(t)o denotes the segment of the path £(t) connecting the points 
E(£, 0) and £(t, s), and W[ s £(t) ] is the corresponding (open) Wilson line. 
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Finally, Pt in eq. ( |1 1 1| ) denotes ordering of the t integral, while the s integral 
is not ordered (although there is an s-ordering inside each M / [ s E(t) ].) 

In the abelian case the two open Wilson lines VT[ s E(t) ] in eq. ( |111| ) cancel 
each other and one recovers the usual Stokes theorem, involving only the 
gauge field strength. In the non-abelian case, however, the factors W[ s E(i) ] 
are needed for gauge invariance, and introduce an explicit dependence of the 
Wilson loop operator on the gauge potential A^. Thus, as opposed to the 
abelian case, the non-abelian Wilson loop operator cannot be mapped in a 
straightforward way to a fermionic operator through the bosonization rule in 
eq. flTU9|). 

For planar loops this difficulty is only apparent. Indeed, consider W[<9E], 
with E contained, say, in the (1, 2) plane. Imposing the A% = gauge 
condition, there is a remnant gauge freedom for the Ai and A 2 components 
in the (1,2) plane, which is the symmetry of a 2-dimensional gauge theory in 
that plane. As discussed in p0[, one can use that gauge symmetry, together 



with the freedom of parametrization of the surface E, so the open Wilson 



line elements in the right hand side of eq. (|1 1 1| ) become the identity. More 
precisely, choosing the gauge condition A 2 = on the E-plane, eq. ( |111| ) can 
be simplified to 

rl rl G Y\ U 

W[dZ] =trP t exp{z/ dt ds- _F^(E(t, S ))} (112) 

Jo Jo OS ot 

provided that E is parametrized so as to have dH/dt and <9E/<9s parallel 
to the Xi and x 2 axis, respectively. This apparent breaking of rotational 
invariance, which includes the presence of t-ordering but not of s-ordering, 
is a consequence of the A 2 = gauge condition on the E-plane, and will be 
removed by the functional integral over the gauge fields. One more ingredient 
is need: by an appropriate addition of a BRST exact form, the bosonization 
rule ( |109| ) can be written in the covariant way 

$V&^ - ±-±-e lwa K x (113) 

Then, writing 

J[E] =trP 4 exp{±167r/ dt ds- _ e A j A [£(t,s)]} (114) 

Jo Jo OS ot 

with in the fermionic current , the bosonization formula gives 

{J[T\) ferm = (W[dT])cs (H5) 
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where in the left hand side the subindex 'ferm' stands for free fermions. 
This is the non-abelian generalization of the result obtained in [|l|. It relates 
a suitably defined non-abelian flux of the fermionic current through a flat 
surface, and the Wilson loop associated to the boundary of that surface, 
with both quantities in the same representation of the group. 

It should be stressed that the bosonic side of this relation is, by definition, 
independent of the surface S and its parametrization. In the fermionic side, 
however, this is not obvious. The relation was derived assuming a flat surface 
E, and it is tempting to assume that this may be extended to smooth defor- 
mations away from the plane. But more relevant is the apparent breaking 
of rotational invariance in the fermionic side due to the remaining t-ordering 
in eq. ( |114| ). This should certainly be expected to be taken care of by the 
particular parametrization assumed above for S. Indeed, one should expect 
that the very need of a parametrization and of a matching ordering of the 
surface integral of the fermionic current, is just a limitation of our present 
analysis. In addition, as is well known, the expectation value of the Wilson 
loop is singular and must be regularized. A natural and consistent regular- 
ization scheme is provided by the framing of the loop f31 |. In the case of a 



non-intersecting loop on a plane, considered here, that framing can be chosen 
also as a plane loop, not intersecting itself nor the original loop. Again, it is 
not clear at this point how these singularities in the bosonic side will show 
up in the (free) fermionic side, and what role will the framing play from the 
fermionic point of view. 

It is natural to ask whether this analysis can be extended to several loops 
and their possible linkings, as done in for the Abelian case. In the bosonic 
side one is interested in the expectation value (W / [ri]W / [r 2 ])c5 or, better yet, 
in the ratio 

, <^™% (U6) 

For non intersecting loops this is a well defined, non singular object in the 
Chern-Simons theory, which depends only on the linking of the two loops Ti 
and T 2 [5I|. Assuming this to be non-trivial (and non-singular), the two loops 



cannot be flat and lying on the same plane, so the previous construction fails. 
But once the ratio ( |116| ) has been computed in the Chern-Simons theory, 
we can take the limit in which the two loops collapse onto a single plane. 
This is a singular limit in which the loops necessarily intersect each other. 
Their linking is not well defined any more, and the value of ( |116| ) depends 
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Figure 1: Different overlaps of the surfaces Ej and E 2 on a plane, determined 
by the possible liftings of the loops <9£i and <9E 2 away from the plane. 



on the initial non-singular loops used in the computation. However, at the 
classical level, before the functional integral is performed, we can repeat the 
previous construction with no difficulties for any arrangement of loops on the 



plane j4(J. Thus, formally we can write 

(J\Ex U E 2 ]) /erm (WidZtjWid^cs 



(J[^l])ferm(J[^2])ferm {W^}) CS (W[dT, 2 ]) CS 



[in) 



where both surfaces Ex and E 2 are contained in the same plane. As we just 
stated, the bosonic side of this relation will be ill defined in general. But it can 
be given a well defined meaning by lifting the loops <9Ej from the plane to non 
intersecting three-dimensional loops Tj. This can be done in different ways, 
specifying different linkings of the loops Tj compatible with the intersections 
of their projections <9Ej onto the plane. Correspondingly, in the fermionic 
side, the surface Si U E 2 must be complemented with a prescription stating 
the way in which the two surfaces overlap. The different possible liftings of 
the loops specify different overlaps of the surfaces, as ilustrated in Fig. (1). In 
this way, relation (|117|) (and its generalizations) can be viewed as a defining 
relation, through bosonization, of the vaccuum expectation value of the flux 
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of the fermionic current through surfaces with foldings. 



8 Bosonization in d = 4 

In ref.|| we have applied our bosonization approach to the study of vector 
and axial-vector currents in d = 4 dimensions. I will briefly describe the 
main results of this work. Consider the generating functional for a massless 
Dirac field in 3 + 1 (Euclidean) dimensions, coupled to Abelian vector (s^) 
and axial- vector (£ M ) external sources 

Zis^t^) = J VipVip exp -S^.^s^t^) 

S$, ^ s M , Q = -i J d 4 x $ - 75^) ^ (118) 

where we have adopted for the 7-matrices the following conventions: 

l\ = 7m , 75 = 75 , {7m , lu} = 2 6^ . (119) 

The addition of the source s M is due to the fact that, in four dimensions, the 
vector and axial currents are independent fermionic bilinears. Thus not all 
the information provided by Z(s^, t^) can be obtained from, say, Z(s fl , 0). In 
two dimensions, because of the smaller number of generators for the Dirac 
algebra, these two currents are related, and the bosonization rule for one of 
the currents also implies the proper rule for the other. 
We now consider the following change of variables 



In terms of the new variables, the generating functional ( |118| ) reads 

(121) 

where the primed fermionic fields have been renamed as unprimed for the sake 
of simplicity, and J is the anomalous Jacobian corresponding to this fermionic 
change of variables, a well-known consequence of the chiral anomaly [f|l 
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This Jacobian is evaluated by using the standard Fujikawa's recipe accomo- 
dated so as to get the consistent anomaly. The answer is [B] 



J (ci, Sp, t^) 



exp 



1 



47T 



■ (122) 



Now, as ( |120| ) is a change of variables, Z cannot depend on either 9 or a. 
Thus 9 and a can be integrated out without other effect than the introduction 
of an irrelevant constant factor in Z, which we ignore. Integration over 9 and 
a is equivalent to integration over two flat Abelian vector fields 9p and a p. 



9,, 



9^9 , ap 



dp a 



f^u{9) = dp9 u - d u 9p = , fp V {a) = d^a u - d v otp 



Writing 

J(ap- 7 Sp,tp) = exp 
eq. (|121|) becomes 



4vr 2 



f 4 1 

J d X Oip^x) tpv pai^S v d pSu -\- ~t u Opt a/ 



(123) 



(124) 



Z = J V9pVapV^V^ 5[fp U (9)}5[fp U (a)} %; S/1 ,i, 



exp 



■S(i/>,ip; Sp + 9p,tp + otp) 



(125) 



Formally integrating out the fermionic fields and making the shift of variables 



6 p 9 p Sp , dp ► (Xp tp , 



(126) 



(IHD leads to 

Z(-Sp,tp) = J V9pVa>p5[fp V (9 - s)}5[fp U (a-t)}J(a>p-tp]Sp,tp) 

x det(^ + i ^ + «7 5 tf) . (127) 

As before, we exponentiate the two functional delta functions in (|127|) using 
two antisymmetric tensor fields Ap V and Bp V as Lagrange multipliers 
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x exp (i J d A x[e^ p(7 A^{d p e (J - d p s a ) + e pvp(T B pu (d p a a - d p t a )} 

x det(^ + z ^ + i7s ^) . (128) 
The bosonized form of Z can then be obtained by integrating out P and a p 



in 



128|) . This produces a generating functional with the tensor fields A 



and B pv as dynamical variables. This step requires the evaluation of the 
fermionic determinant, which of course is necessarily non-exact. 

At this stage we can already derive the rules that map the vector and 
axial- vector currents into functions of the bosonic fields A pu and B pv . This 
correspondence requires no approximation and may well be called 'exact'. 
These rules follow from elementary functional differentiation 

3 ix = (^7^) = ^J^T logZ| SM=0 = -e^ pa d„A pa (129) 

5 _ i 
kn = (^757^^) = -«^-log^| tM=0 = -e pupa d u B pa - — e pupa s u d p s a . 

(130) 

From the antisymmetry of the tensors A pv and B pu , we are entitled to derive 
the equations for the divergencies of the currents: 



d»3n = 



i 



W = -^ F M F M- ( 131 ) 



with F pv = (l/2)e pva pF a p. We then see that the bosonization rule ( |130| ) 
correctly reproduces the axial anomaly. 

As before, although the bosonization recipe ( |129| )-( |T30|) for associating 
the fermionic currents with expressions written in terms of bosonic fields 
is exact, the bosonic action governing the boson field dynamics cannot be 
evaluated in an exact form in d > 2 dimensions. Different approximations 
for computing the fermionic determinant would yield alternative effective 
bosonic actions valid in different regimes. In ref . || the fermionic determinant 
in ( p.28| ) was evaluated to second order in the fields P and a p . The use of this 
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quadratic approximation can be motivated by the same kind of arguments 
(see in particular the 'quasi-theorem') used in ref. |1J. Calling 

det(0 + %9 + ij 5 $ = exp [W(9^, eg] . (132) 

the answer for the renormalized W is 

1 



W{6^ap) = -- I d 4 xd A y 9^x)6^F(x - y)9 v {y) 



+ <*p(x)6±,G(x - y)a v (y) + m 2 a^x)^ 5(x - y)a u (y) 



where 



G{x — y) = F(x — y) + m 2 5(x — y) . 
With this result, one easily finds for the generating functional 



(133) 
(134) 



Z(s fl ,t fl ) = exp[C(s M ,* M )] f VA^VB^ x 

exp ^ i J d x\s p j €p J i,p a dpAp (7 -\- 1 pitpvpcfdpB p a -\- ~^^tp U p (T SpdpS a y\^ x 
expj-i J d 4 xd i y[Ap Up (x)F~ l (x - y)A fa/p (y) + 
Bp Vp (x)G~ l (x - y)Bp Vp (y)}} x exp j-^ J d 4 xd 4 yd IM B up (x)x 



G 1 (x - y)5 fll/p>a/37 (s a dps 1 + -tadpt-,) 



(135) 



where 



A 


— dpA up + d u A P p 


+ dpApj, 


B p,vp 


= dpB vp + dyBpp 


+ dpBpy 






' Spa 


<Vr \ 


pup,a/3-y 


= det 


Sua u p 








\ Sp a 5pp 


5 P1 ) 



(136) 



and 



2(2tt) 



d xd y \ [sp(x)d l/ s x (x) + -tp(x)d l/ t x (x)} 
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J d 4 xd 4 yG(x)d- 2 G-\x-y)g(y) 



2(2irY 

+ 2m2 [ 2?r)4 / #xd*yG(x)dr*(x -y)Q{y)} (137) 

where Q = e^ p \(d^s„d p s x + \d p t v d p t x ). 

In conclusion, we have applied our bosonization technique to the case of 
massless Dirac fermions in four dimensions in the presence of both vector and 
axial-vector sources. This has allowed us to find the bosonization rules for 
both fermionic currents, eqs.( |129[ )-( |13U| ), in terms of Kalb-Ramond bosonic 
fields. While the bosonization rule for the vector current can be written in 
a natural and compact form, reminiscent of the well-known two-dimensional 
bosonization rule, 

^li^ -£^ P ad u A pa , (138) 
the result for the axial current is more involved and includes the vector source 

_ 2 

V>757/^ -> -ttxvpadvBpv - — tyvpaSydpSa . (139) 

In our approach, this is a consequence of the anomalous behaviour of the 
fermionic measure under axial gauge transformations and in this way the 
bosonic form of the axial current correctly yields its anomalous divergence. 
We want to comment on the possibility of considering the particular case of a 
purely chiral external source (s^ = ±t^), and obtaining a bosonized version 
for this model. The Kalb-Ramond field then corresponds to a particular 
'chiral' combination of A and B, namely A^ u ± B pu . 
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